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ABSTRACT 

We have conducted an archival Spitzer study of 38 early-type galaxies in order to determine the origin of the 
dust in approximately half of this population. Our sample galaxies generally have good wavelength coverage 
from 3.6/im to 160/im, as well as visible- wavelength HST images. We use the Spitzer data to estimate dust 
masses, or establish upper limits, and find that all of the early-type galaxies with dust lanes in the HST data 
are detected in all of the Spitzer bands and have dust masses of ~ 1O 5_6 - 5 M , while galaxies without dust lanes 
are not detected at 70/im and 160/im and typically have < 1O 5 M of dust. The apparently dust- free galaxies 
do have 24/im emission that scales with the shorter wavelength flux, yet substantially exceeds the expectations 
of photospheric emission by approximately a factor of three. We conclude this emission is dominated by 
hot, circumstellar dust around evolved stars that does not survive to form a substantial interstellar component. 
The order of magnitude variations in dust masses between galaxies with similar stellar populations rules out 
a subtantial contribution from continual, internal production in spite of the clear evidence for circumstellar 
dust. We demonstrate that the interstellar dust is not due to purely external accretion, unless the product of the 
merger rate of dusty satellites and the dust lifetime is at least an order of magnitude higher than expected. We 
propose that dust in early-type galaxies is seeded by external accretion, yet the accreted dust is maintained by 
continued growth in externally-accreted cold gas beyond the nominal lifetime of individual grains. The several 
Gyr depletion time of the cold gas is long enough to reconcile the fraction of dusty early-type galaxies with the 
merger rate of gas-rich satellites. As the majority of dusty early-type galaxies are also low-luminosity Active 
Galactic Nuclei and likely fueled by this cold gas, their lifetime should similarly be several Gyr. 
Subject headings: dust - galaxies: general - galaxies: ISM - infrared: galaxies - ISM: general 



1. INTRODUCTION 

One of the most striking characteristics of elliptical and 
lenticular galaxies is their apparent uniformity. Much more 
than their later- type cousins, early-type galaxies appear to 
have quite smooth and symmetric surface brightness profiles. 
Their morphological self- similarity is largely because most 
early-type galaxies, and in particular ellipticals, are expected 
to be the result of nearly equal-mass mergers that have pro- 
duced symmetric stellar distributions via dynamical relax- 
ation. This assembly process, combined with the generally 
large halos that host early-type galaxies, has also substan- 
tially heated the gas in their interstellar medium (ISM) such 
that cooling and subsequent star formation are quite ineffi- 
cient. As a result, their stellar populations tend to be mostly 
old and well-mixed, which reinforces their similar and smooth 
appearance. 

Against the backdrop of these similarities, and the simplic- 
ity of this general picture, the differences between nominally 
similar early-type galaxies are cast into sharper relief than for 
the more highly structured late-type disk galaxy population. 
One striking difference is that some, but not all, early-type 
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galaxies clearly contain cold, interstellar dust. This was re- 
vealed in IRAS observat ions that detected a number of ellip- 
tical galaxies at lOQmn (jJura et al.l[l9 87: Knapp et alJfl989b 
Goudfro oii & de Jong] 1 19951: iBregman et al. 1998). A study 
by iKnapp et all d!992|) also demonstrated that ellipticals had 
emission at 10- 12/nm in excess of what was expected from 
stellar photospheres alone. This excess could be explained 
by circumstellar dust emission from evolved stars and proved 
to be a strong observational confirmation of expectations for 
mass loss from old stellar populations. 

The presence of dust in early-type galaxies is interesting 
because it may provide additional information about their for- 
mation histories, the properties of dust when embedded within 
a substantially greater fraction of hot gas than is present in 
typical spirals, and the evolution of the interstellar medium, 
in addition to stellar mass loss from evolved stars. Yet the 
presence of dust is also in apparent conflict with the expec- 
tation that the dust destruction timescale should be relatively 
short (jdust ~ 2 x 10 4 y r), largely due to sputteri ng of small 
grains in their hot ISM (iDraine & Salpetedll979bl) . 

Observations with HST revealed that dust lanes were com- 
mon and probably present in the majority of early-type 
galaxies. These observations showed dust lanes in ab- 
sorption within 100s of parsecs to a kpc of the nucleus 
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( van Dokkum & Franx 1995). This study also noted that the 
incidence of dust is higher is radio galaxies than in radio-quiet 
galaxies, estimated that the mass of dust ranged from 10 3-7 
M , and that the distribution of the dust lanes sugg ests exter- 
nal ac cretion. Subsequent studies of atomic (Morganti et al] 
120061) and mole cular gas kinematics (I Young etaLl 120111: 
iDavis et all 1201 1[) supported external accretion because the 
gas kinematics are often decoupled from the stellar kinemat- 
ics. 

While there is good, anecdotal evidence for external accre- 
tion, the evolved stars in e arly-type galaxies ar e also a po- 
tential source for the dust dKnapp et al.lll992|) . lAthey et all 
( 2002) showed that the mid-infrared emission from early-type 
galaxies is consistent with stellar mass loss of on order 0. 1 - 1 
M yr -1 for a typical L* galaxy. For a dust-to-gas ratio of 
e dg = 0.005 and a dust lifetime of 10 7-8 yr characteristic of the 
Milky Way, this would lead to a typical reservoir of 5 x 10 3-5 
M (or 10-100 M for r dust = 2 x 10 4 yr). While this dust 
would initially match the spatial and kinematic distribution of 
the evolved stellar population, Mathews & Brighenti (2003) 
suggest that cooling in dust-rich gas could lead to clumping 
and settling into the central kiloparsec on a timescale of 10 8-9 
yr. While this substantially exceeds the expected lifetime of 
individual dust grains, the dumpiness may sufficiently shield 
the dust to produce a steady state dust mass that is consistent 
with the observations. However, the internal origin hypothesis 
does not explain the dramatic vari ation in dus t mas s between 
otherwise quite similar galaxies. iTemi et al] (|2004|) showed 
that there is no correlation between the far-infrared luminos- 
ity of early-type galaxies and their visible-wavelength lumi- 
nosities, even though visible-wavelength light should corre- 
late reasonably well with the number of mass-losing stars. 
This absence of a correlation can be more readily explained 
by the external accretion hypothesis. 

Another important aspect of these studies was the obser- 
vation that the majority of the early-type galaxies with dust 
also have emission-line nuclei that o ften indicate the pres- 
ence of a LINER or Seyfert galaxy (|va n Dokku m~& Franxl 
119951: iRavindranath et alJl2001t Lauer et al.l 12005b . In order 
to investigate this corre lation with a well- selected sample of 
early-type galaxies, in ISimoes Lopes et all (|2007|) we stud- 
ied archival HST images of a carefully-matched sample of 
active and inactive early-type hosts and determined that all 
active early-type galaxies also have dust within 100s of par- 
sees of their nuclei, while only ~ 25% of inactive galax- 
ies had evidence for dust. This study showed that ~ 60% 
of the early-ty pe galaxy population has interstellar dust. As 
we noted in ISimoes Lopes et all (l2007|) . the high incidence 
of dust in early-type galaxies presents a key challenge to the 
external origin hypothesis, namely that ~ 60% of the early- 
type galaxy population must have accreted a gas and dust-rich 
dwarf within a time period comparable to the dust destruction 
time. For a dust destruction timescale of 10 7-8 yr, this implies 
a very high rate of gas-rich mergers. 

A missing ingredient in the ISimoes Lopes et all ( 2007) 
study is that there was no constraint on the dust mass distribu- 
tion in the early-type galaxies, which would provide a useful 
constraint on the mass range and consequently mass ratio of 
the satelli te galaxies required to p rovide the observed mass 
of dust. Simo es Lopes e t al. (20 071) also emp loyed the struc- 
ture map technique of lPogge & Martini] (|2002|) to identify dust 
lanes. As this is a contrast enhancement technique, it is in- 
sensitive to uniform or diffusely distributed dust. Uniformly- 



distributed dust in the galaxies without dust lanes in the cen- 
tral kpc would provide strong support for the internal origin 
hypothesis. 

Over the last few years, several studies with Spitzer and 
Herschel have assembled far-infrared measurements and esti- 
mated dust masses and basic dust properties for many early- 
type galaxies. One result of this work is that the dust in early- 
type galaxies is warmer than in spiral galaxies dBendo et al.l 
[20131: ISkibba"etani20TTl: ISmith et al.ll2012l: lAuld et al.ll2013l) . 
which may be due to more intense radiation fields or differ- 
ent dust grain properties. The Herschel KINGFISH survey 
showed that early-type galaxies exhibit a strong correlation 
between the dust to stellar flux ratio and the specific star for- 
mation rat e, which could be d ue to low levels of on-going star 
formation ( Ski bba et al.ll201 1|) . The Herschel Reference Sur- 
vey ([Smith et al] 120 1 2l) and the Herschel Virgo Cluster Sur- 
vey (lAuld et al.l 1201 3|) detect many early-type galaxies and 
measure du s t mas ses in the range of ~ 10 5-7 M . The 
Smit h et aT] ([20 1 2 ) study concluded that much of this dust 
was acquired from interactions due to the wide range in dust 
to stellar mass ratio. A similar conclusion was reached by 
Rowland s et al] (l2012|) . who detected 5.5% of luminous ETGs 
in the Herschel AILAS/GAMA study. While only sensitive 
to quite high dust masses (their mean detected dust mass was 
5.5 x 10 7 M ), they conclude that these high dust masses 
could not be produced internally. 

In this paper we present a Spitzer archival study of the ori- 
gin of the dust in earl y-type galaxies. Our samp le is com- 
posed of galaxies in the Sim oes Lopes et al.l (|2007|) early-type 
galaxy sample that have both IRAC and MIPS observations. 
We describe this sample in §|2] and our data processing and 
photometry in §|3] All of the galaxies in this sample are de- 
tected in the IRAC bands, which are dominated by stellar 
emission, as well as in the 24/im MIPS band. Many are also 
detected in the longer wavelength MIPS bands. The 24/im 
emission from galaxies not detected at longer wavelengths 
appears to be dominated by circumstellar dust. We describe 
how we account for this emission in §@] To derive the dust 
masses and proper ties of the dust, we fit these galaxies with 
the dust models of Dra ine & Lil (|2007|) . These fits and the de- 
rived properties are described in §|5] We use these data to de- 
velop a new hypothesis for the origin of the dust in early-type 
galaxies in §0 and summarize our results in §|7] 

2. SAMPLE 

Our sample is composed of 38 bright, nearby, early-type 
galaxies (types E, E/S0, and SO). They constitute all of the 
early- type galaxies from the HST study of Sim oes Lopes et al.l 
( 2007) that have observations in at lea st two MIPS bands. The 
goal of the ISimoes Lopes et al] (l2007[) study was to compare 
dust structure in the central 100s of parsecs between active 
and inactive galaxies and that sample contained equal num- 
bers of well-match ed active and ina ctive galaxies from the the 
Palomar survey of HoIetaD (dUl. The Palomar survey is 
very well suited for this purpose because it consists of very 
high signal-to-noise ratio (SNR) spectroscopy of all bright 
galaxies in a fixed region of the sky. Furthermore, all of the 
nuclear spectra have been classified as either absorption-line 
nuclei, H II galaxies, LIN ERs, transition objects, or Seyfert 
galaxies dHo et al.lll997b|) . The Palomar survey sample is 
thus both unbiased with respect to nuclear activity and has 
homogeneous and sensitive nuclear classifications. A further 
advantage is that there are a lar ge numbe r of ancillary data 
products, which are described in Ho (2008|). 
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TABLE 1 
Sample Properties 



Fig. 1 . — Sample distribution as a function of morphological type, activity, 
and distance. (Top panel) Active and inactive galaxies (red and blue cir- 
cles, respectively) as a function of morphological type and distance. (Bottom 
panel) Histogram of all morphological types of the active and inactive galax- 
ies (red and blue lines, respectively) as a function of distance. The substantial 
overdensity at 16.8 Mpc is because many galaxies in that bin are members of 
the Virgo cluster. 



The ISimoes Lopes et all (|2007|) study began with all of 
the LINERs and Seyfert galaxies (hereafter simply active 
galaxies) in the Palomar sample with HST observations and, 
for each active galaxy, identified a control (inactive) galaxy 
with the same morphology, distance, luminosity, and axis ra- 
tio. They identified 26 early-type pairs with relatively strict 
matching criteria (morphological T type within 1 unit, ab- 
solute B magnitude within 1 mag, distance within 50%) and 
eight additional pairs of early-type galaxies that were not as 
well matched. These 34 pairs of active and inactive early-type 
galaxies had very similar distributions of distances, luminosi- 
ties, and morphological types and are consequently very well 
suited to compare the differences between active and inac- 
tive early-type galaxies. As noted in the Introduction, that 
HST study found that all active, early-type galaxies possess 
some form of dust structure within several 100 pc of their 
nuclei, whereas only 26% (nine) of the inactive galaxies ex- 
hibited evidence for such dust. These data thus show there 
is a strong dichotomy between dusty and active early-type 
galaxies and typically dust- free and inact i ve early -type inac- 
tive galaxies. Furthermore, as iHo et al.l (Il997b|> found that 
50% of all early-type galaxies are classified as LINERs or 
Seyferts (active), this result also indicates that about 60% of 
all early-type galaxies have dust. 

We construc ted our sample from a ll of the early-type galax- 
ies in the Simoes Lope s et all (120071) study that have observa- 
tions in at least two MIPS bands, and in most cases all three 
are available. All of the galaxies also have IRAC observa- 
tions, in addition to the HST observations. This sample has 
approximately equal numbers of active a nd inactive galaxies 
and sh ould be an unbiased subset of the Sim oes Lopes et al.l 
(|2007|) sample because these galaxies were targeted for MIPS 
observations by a multitude of small programs. Because of 
this selection, and a desire to not decrease the sample size 
further, we did not attempt to create a pair-matched sample 
from this subset. We did confirm that there are not obvious 
trends between the active and inactive subsets as a function 
of distance and morphological type. Figure Q] shows the dis- 
tance distribution for galaxies classified as type E, E/S0, and 
SO. While the true ellipticals extend out to larger distances 
than the two other types, the distribution of active and inactive 
galaxies are similar. The morphological classification, activ- 



Name 


Morph 


Activity 


Distance 


HST Dust? 


(1) 


(2) 


(3) 


(4) 


(5) 


JNCjCU315 


E 


active 


65.8 


Y 


NGC0821 


E/S0 


inactive 


23.2 


N 


NGC1023 


SO 


inactive 


10.5 


N 


JNCjC23UU 


E/S0 


inactive 


31.0 


N 


NGC2768 


E/S0 


active 


23.7 


Y 


JNCjC2/o/ 


SO 


active 


13.0 


Y 


JNCjC322o 


E 


active 


23.4 


Y 


JNCjC33 / / 


c /o a 


inactive 


8.1 


Y 


JNCjC3414 


SO 


active 


24.9 


Y 


JNCjC36U7 


SO 


active 


19.9 


Y 


NGC3640 


E 


inactive 


24.2 


N 


NGC3945 


SO 


active 


22.5 


Y 


NCjC399o 


so 


active 


21.6 


Y 


JNCrC4UZo 


OA 


inactive 


1 1 A 
1 /.U 


Y 


JNCjC413o 


OA 


active 


1 1 A 
1 /.U 


Y 


JNCjC42/o 


E 


active 


9.7 


Y 


NGC4291 


E 


inactive 


29.4 


N 


XT/" 1 /" 1 A 

NGC4293 


SOa 


active 


17.0 


Y 


"XT/" 1 /" 1 A 1£LZ 

JNCjC43o5 


E 


inactive 


16.8 


N 


"VT/" 1 /" 1 A T7 1 

JNCjC4371 


SO 


inactive 


16.8 


Y 


XT/" 1 /~i/(TOO 

NGC4382 


so 


inactive 


16.8 


N 


XT/" 1 /" 1 /! Af\/1 

NGC4406 


E 


inactive 


16.8 


N 


"XT/" 1 /" 1 A CO/; 

JNCjC452o 


SOa 


inactive 


16.8 


Y 


M/" 1 /" 1 /! CCA 

JNCjC4jj(J 


SO 


active 


16.8 


Y 


M/" 1 /" 1 /! C7A 

JNCjC4j /(J 


SO 


inactive 


16.8 


N 


"VT/" 1 /" 1 /! CTO 

NGC4578 


SO 


inactive 


16.8 


N 


NGC4589 


E 


active 


30.0 


Y 


NGC4612 


SO 


inactive 


16.8 


N 


NGC4621 


E 


inactive 


16.8 


N 


NGC4636 


E 


active 


17.0 


Y 


NGC4649 


E 


inactive 


16.8 


N 


NGC4694 


E 


inactive 


16.8 


Y 


NGC5077 


E 


active 


40.6 


Y 


NGC5273 


SO 


active 


21.3 


Y 


NGC5308 


SO 


inactive 


32.4 


N 


NGC5557 


E 


inactive 


42.6 


N 


NGC5576 


E 


inactive 


26.4 


N 


NGC7619 


E 


inactive 


50.7 


N 


Note. — Properties of the 


galaxies in the sample. 


Columns are 


: (1) Galaxy 



name; (2) Morphological type; (3) Activity type; (4) Distance in Mpc; (5) Presence 
or absence of dust in HST data from Simoes Lopes etal. (2007) References for 
morphology, activity, and distance are provided in Simoes Lopes et al. ( 2007). 

ity type, and distance for all of the galaxies are presented in 
Table[U along with whether o r not they were observed to have 
dust structures in the study of Simoes Lop es et ail (|2007). 

3. DATA PROCESSING AND PHOTOMETRY 

We obtained all available IRAC and MIPS data for these 
galaxies from the Spitzer archive. The IRAC data correspond 
to images in bandpasses centered at approximately 3.6/im, 
4.5/im, 5.7/im, and 8/im. For most galaxies, the first three 
bands are dominated by stellar emission, while the fourth is 
dominated by emission from polycyclic aromatic hydrocar- 
bons (PAHs). For most early-type galaxies, the fourth band 
is also dominated by stellar emission, rather than PAHs, as 
we discuss further below. The MIPS data correspond to im- 
ages in bandpasses with nominal centers at 24/im, 70/im, and 
160/im. Emission at 24/im is typically dominated by hot dust, 
which is either interstellar or circumstellar, while emission at 
longer wavelengths is dominated by cooler, diffuse, interste- 
lar dust. We typically started our analysis with the Basic Cal- 
ibrated Data (BCD) for the MIPS observations and the Post- 
BCD (PBCD) products for the IRAC observations. In the sub- 
sections below we describe any additional processing, photo- 
metric measurements, and how we determined upper limits. 
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Fig. 2. — Structure maps (left), PAH images (middle), and MIPS 160/im images (right) for the sample. The structure map is 18" on a side, while the PAH and 
MIPS images are 180" on a side. The small box on the PAH and MIPS images indicates the field of view of the structure map. The black, dashed circle on the 
PAH image indicates the photometric aperture used for the 3.6/mi to 70/xm photometry (a PSF model was used for the 160/xm data). The third panel is empty if 
160/im MIPS data were unavailable for that galaxy. 



3.1. IRACData 

The PBCD data for these galaxies are sufficiently robust 
that no additional processing was necessary in most cases. 
The few exceptions were galaxies that had relatively few 
IRAC observations. In these cases we obtained the BCD data 
from the archive and processed them into mosaics with the 
MOPEX software package from the Spitzer Science Center. 

We measured aperture photometry in all four IRAC bands 
with the phot task in IRAK We typically adopted an aper- 
ture radius of V because this includes nearly all of the ob- 
served dust emission at longer wavelengths. In some cases 



we used an 0.5 ' radius due to either nearby, bright objects or 
the extent of the imaging data in one or more bands. In sev- 
eral cases bright, foreground stars were in the aperture and we 
masked them out by replacing those pixel values with the me- 
dian of nearby pixels. The aperture size and photometry for 
each galaxy is listed in Table [2 These measurements include 
the extended source aperture correction. 

The apertures are smaller than the full extent of the visible 
emission of these galaxies (e.g. D25). While these smaller 
apertures have less noise, in some cases they may lead to un- 
derestimates of the total emission. 
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NGC4278 



C4293 




NGC4365 



NGC4371 



NGC4382 







160^m 



160/im 



160^171 



160/im 




.0 M m 



Fig. 3. — Continuation of Figure |2] 

3.2. MIPS Data 

The Spitzer MIPS data were processed further with the 
MOPEX package. In addition to the mosaicing and other pro- 
cessing (e.g. median filtering), we used the IDL program BCD 
col umn filter. pro to calculate and subtract the median 
value from each column for each BCD. This process signifi- 
cantly diminishes artifacts in the images with minimal loss of 
flux from the source. 

The fluxes at 24/im and 70/mi were measured with aperture 
photometry in the GAIA package and the same aperture sizes 
employed for the IRAC measurements. Aperture corrections 
were applied to these measurements. These corrections were 
derived from the point spread function (PSF) of bright point 



NGC4526 



NGC4570 



NGC4589 



NGC4612 



160^m 



16(Vm 




160^m 



sources in MIPS data. Aperture photometry proved impos- 
sible to use for the 160/im fluxes because it was difficult to 
accurately estimate the background. This is because the PSF 
at 160/im is comparable in size to the field of view in many 
cases. There is also the strong possibility of contamination by 
other sources. We therefore extracted the 160/im fluxes with 
the APEX package in MOPEX, which uses a PSF-fitting tech- 
nique. The standard PSF available with the MOPEX package 
was used for this task. 

We quantified several sources of measurement uncertainty. 
First, we calculated uncertainties for each galaxy in each band 
from the standard deviation of six measurements of the back- 
ground. For some of the 70/im data and most of the 160/im 
data, this technique may not accurately capture the uncertain- 
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NGC4636 



NGC4694 



i 



NGC5273 



NGC5308 



■ b 



160/^m 



9 




160^m 



9 



Fig. 4. — Continuation of Figure |2 

ties as the field of view is too small. In these cases we set 
a lower limit on the uncertainty based on the noise in other, 
similar datasets w ith a larger field of vie w, such as work by 
ITemi et al.1 (120091) and lBendo e7ail (120121). as well as the con- 
fusion limit measurement by iFrayer et all (|2009). Finally, 
we included the instrumental flux calibration uncertainties (of 
4%, 10%, and 12% for 24/zm, 70/im, and 160/im, respec- 
tively). Calibration uncertainties dominate the uncertainties 
for the majority of the detections. 

All of the galaxies were easily detected at 24/im and the 
apertures were centered on this emission. Detections at 70/im 
and 160/im were consistent with the same centroid as the 
24/im emission. If there was not an obvious detection in one 
of these two bands, then an upper limit was measured at the 



NGC5576 



r 




160miti 



coordinates of the 24/im emission. In all cases we quote So- 
upper limits. Because of the size of our aperture, as well as 
PSF-fitting at 160/im, these upper limits are only valid if the 
dust is centrally concentrated. 

Many of these archiv al datasets were in clud ed in previ- 
ous st udies, in particular ITemi et aT] (|2009|) and iBendo et al.l 
(|2012|) . Our fluxes g eneral ly agree to within the uncertain- 
ties with ITemi et ail ([2009), who employed a similar set of 
data processing and calib ration techn i ques. The main differ- 
ences are likely because ITemi et ail (|2009|) employed aper- 
tures that were comparable in size to the visible extent of the 
galaxy. Our apertures are smaller, often by a factor of sev- 
eral, while at 160/im we employed a PSF-based measurement 
rather than aperture photometry. We have larger differences 
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TABLE 2 
Spitzer Photometry 



Name 


Aperture 


3.55/im 


4.49^m 


5.73/im 


7.87^m 


23.68/im 


71.42/im 


155.9^m 


JNCjCU315 


r 


88.0 


55.0 


39.6 


30.8 


no o (a c\\ 


QOA Q /OO A\ 

32U.3 (33.4) 


475 (58) 


JNCjCUozI 


0.5' 


100.9 


60.6 


40.2 


25.1 


o o (r\ 
0.0 ((J. /) 


< 25.9 


< 87.4 


NCjC1023 


r 


330.0 


197.8 


135.9 


82.7 


CO 1 /<"> 1 \ 

52.3 (2.1) 


< 31.5 


< 29.7 


JNCjC23UU 


r 


122.6 


73.0 


49.4 


30.2 


19 (0.8) 


< 16.9 


< 66.6 


JNCjC27o5 


r 


136.1 


83.8 


58.4 


39.6 


32.4 (1.3) 


rill /C 1 c\ 

513.2 (51.5) 


486 (61) 


MPPn O H 

JNCjC2/57 
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Note. — Photometry of the galaxies in the sample. Columns are: (1) Galaxy name; (2) Aperture radius in arcmin for the 3.6/xm 
to 70/xm measurements (a PSF model was used for the 160 /im photometry, as well as the 70/im photometry of NGC3226, which 
is marked with a p superscript); (3-6) Flux in IRAC channels 1-4 in mJy; (7-9) Flux in MIPS channels 1-3 in mJy. Uncertainties 
include statistical and calibration uncertainties. All upper limits are 3<r upper limits. While NGC4636 is clearly detected at 160 /im, 
we were unable to obtain a good photometric measurement. 

with the Ben do et all £2012) measurements and the general 
sense of these differen ces is that our flux m easurements tend 
to be smaller Because | Bendo et all (|2012) used similar aper- 
ture sizes to Temi et al. (2009), the difference is likely be- 
cause [Bendoetai] (120121) based their pipeline on an older set 
of processing and calibration tools. Two of the most signif- 
icant differences we have traced to either the presence of a 
nearby companion (for NGC3226 ) or clumps of emis sion at 
large radii (for NGC4406, see also lGomez e"t~aD2010l) . 

As our apertures are smaller than the full extent of the stel- 
lar emission (although not necessarily the dust emission), our 
measurements will underestimate the flux if the emission is 
more extended. In most cases these differences appear to be 
at most on the order of a few tens of percent. The only signifi- 
cant difference is NGC4526, where our 160jimi measure ment 
is only 60% of the value reported by Ben do et ail (|2012), and 
about a fac t or of t wo below the PACS 160/im measurement of 
lAuld et ail (|2013|) . which has a much larger field of view and 
consequently better background measurement. We discuss the 
im plicatio ns of our measurements on the inferred dust masses 
in 35X21 

4. STELLAR AND CIRCUMSTELLAR EMISSION 



Most of the dust in galaxies is in the interstellar medium. 
This dust is likely composed of a mix of amorphous silicate 
and graphitic grains, as well as small PAH particles. Both 
stars and PAH particles contribute emission to the 8/im IRAC 
band, so it is important to accurately determine the stellar con- 
tribution to the 8/im band in order to measure the PAH com- 
ponent. Similarly, both hot, circumstellar dust around evolved 
stars and hot, interstellar dust contribute emission to the 24/im 
MIPS band, so it is important to determine the circumstellar 
contribution to the 24/im band in order to quantify the warm 
component of the diffuse ISM. In the next two subsections 
we describe how we used our sample of galaxies with no evi- 
dence for interstellar dust to obtain an empirical measurement 
of the stellar and circumstellar emission in the 8/im and 24/im 
bands. 

4.1. Stellar Contribution 

IRAC 8/im observations of most galaxies include contribu- 
tions from both stellar photospheres and PAH molecules. In 
order to study the total flux and spatial distribution of PAHs, 
the conventional practice is to remove the stellar contribu- 
tion by subtraction of a scaled version of the 3.6/im flux, 
where this scale factor is calculated from models of stellar 
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Fig. 5. — Ratio of 8^m to 3.6/im flux (top panel) and 8/im to 4.5/im flux 
(bottom panel) as a function of 4.5/mi luminosity for galaxies with the small- 
est value of this ratio. Red points correspond to galaxies without evidence for 
dust, both based on upper limits at 160/im and the absence of dust structure 
in HST imaging. Blue points correspond to galaxies with 160/im detections 
and dust structure in HST imaging. The galaxies with no evidence for dust 
have median ratios of / 8 // 4 .5 = 0.418 ± 0.011 and h/h.e = 0.251 ± 0.010 
and we adopt these values to remove the stellar contribution at 8/xm for our 
early-type galaxy sample. The errorbar in each panel represents the average 
uncertainty in the measurement of the ratio. Several galaxies with detections 
at 160/im (blue points) have similar value of this ratio, although some are 
substantially higher and off the upper end of the range shown on this figure. 



photospheres. Hel ou et ail (120041) estimated that this value is 
0.232 based on stella r population models from Starburst 99 
(Leit herer et al.lll99"9|) and noted this factor is not very sensi- 
tive to star formation history or metallicity. They do note that 
the 3.6/im emis sion may include a hot dust componen t, based 
on earlier work (iBernard et al.ll994i iHunt et al.ll2002l) but es- 
timate that this is a small contribution. iDraine et al.l (|2007|) 
obtained a similar value of 0.260 with the assumption that the 
emission is described by a 5000K blackbody. 

The stellar contribution at 8/im is a small fraction of the 
total emission from the late- type spiral NGC 300 studied by 
iHelou et ail (H004). This is also true for most of the galaxies 
in the SINGS sample studied bv lDraine etaLl (12007b . While 
approximately half of the galaxies in our sample are similarly 
expected to have substantial dust emission at 8/im, others may 
have little to no dust emission in this band. For these galaxies, 
an accurate estimate of the scale factor specific to early-type 
galaxy stellar populations is necessary to correctly identify 
and measure PAH emission at 8/im. 

We have empirically estimated the appropriate scale fac- 
tor for early-type galaxies with the nominally dust- free galax- 
ies in our sample. Figure shows the ratio of /s//3.6 and 
hi U.5 for the subset of the sample with the smallest val- 
ues of these ratios as a function of 4.5/im luminosity. In 
both cases all of the dust-free galaxies have similar values 
of these ratios, which correspond to fs/h.6 = 0-25 ± 0.01 and 
h/ U.5 = 0-42 ± 0.01. The uncertainty quoted for both ratios 
represents the rms variation of the dust- free galaxies. The 
average photometric uncertainty for the ratio measurement is 
also shown in each panel, although because it includes (corre- 
lated) absolute calibration uncertainties it is an overestimate 
of the true uncertainty in the individual points. The measured 
value for A/jjfi * s con s istent at la with the value of 0.260 
chosen by IDraine et"aI1 (l2007|) and is 2a from the value of 
0.232 calculated bv (IHelou et ai1l2004l) . 

We have compared these values to the predic t ed flu x ra- 
tios from stel l ar po pulation models by iMarastonl (120051) and 
IConroy et al.l (|2009l) to judge how well the models fit these 



flux ratios, examine the expected variation due to stellar pop- 
ulation differences, and estimate the significan ce o f the 24/im 
detections of apparently dust- free galaxies (in §14. 2t . As early- 
type galaxies are primarily composed of old stellar popula- 
tion, and this component is expected to dominated the stellar 
emission in the infrared, we only examine models with single 
stellar population a ges, but consider metallicites from solar to 
super solar. From Marastonl (|2005|) we calculate the fo/fas 
and fs/fa.s ratios for metallicities of Z = 0.02, Z = 0.04, and 
Z = 0.07, a Kroupa IMF, and t h eir red horizontal branch pre- 
scription. From IConroy et al.l (|2009) we calculate these ra- 
tios for Z = 0.019 and Z = 0.035, a Kroupa IMF, the Padova 
isochrones, and the Basel library. 

The flux ratios for these models as a function of popula- 
tion age are shown in Figure |6] The fs/h.6 and fs/f 4.5 ratios 
from the model s agree we l l with the empirical data shown in 
Figure [5] The IMarastonl (|2005|) models indic ate a trend of 
larger fs/ h.6 ratio at larger metallicities. The IConroy et al.l 
( 2009) models appear relatively insensitive to metallicity, but 
do show some increase at larger ages. The fs/f 4.5 ratios for 
both sets of models agree well with each other and the data, 
and also appear essentially independent of age and metallic- 
ity. Models with different IMF prescriptions lead to essen- 
tially identical results and are not shown. 

4.2. Circumstellar Dust 

While the galaxies without evidence for dust in HST im- 
ages were also undetected at 70/im and 160/im, all are de- 
tected at 24/im. Figure display s a similar plot to Figure [5] 
but for 24/im rather than 8/im. This figure indicates that the 
24/im emission is relatively constant for the apparently dust- 
free galaxies, although the scatter is larger than for the 8 /im 
emission. The larger scatter may be due to stellar popula- 
tion differences, such as age or metalli city. This was indi- 
cated by the work of lAthey et al.l (|2002|) . who measured sig- 
nificant intrinsic scatter between the mid-infrared and B-band 
luminosities of nine early-type galaxies. The ratios and the 
rms variation we measure are /24/A5 = 0.253 ±0.074 and 
/W/3.6 = 0.154 ±0.047. While some of the dusty galaxies 
also have consistent values of these ratios, the ratios for many 
are much larger and some are not shown on the figure. 

These ratios for all galaxies, including the apparently dust- 
free galaxies, are several times larger than expected from stel- 
lar photospheric emission alone. The right panels of Figure [6] 
show the empirical .fa/.h .6 an d fa I A.5 ratios com pared to 
the same IMarastonl (120051) and IConrov et all (120091) models 
described in the previous section. In both cases the empiri- 
cal values are approximately a factor of three higher than all 
of the models. The models are also in good agreement with 
each other and show little evidence for significant variation as 
a function of stellar population age or metallicity. We there- 
fore conclude that the 24/im emission is clearly in excess of 
stellar photospheric emission, yet it does scale with the flux 
of the stellar population. 

The most likely origin of this excess 24/im flux is emission 
from hot dust in the circumstellar envelopes of evolved stars. 
For a population of galaxies with relatively uniform and old 
stellar populations, the fraction of evolved stars with circum- 
stellar dust should be an approximately fixed fraction of the 
stellar mass, which in turn is reasonably well represented by 
the flux at 3.6/im and 4.5/im. Mod els of infrared emi ssion 
from dust shells around AGB stars by iPiovan et al.l (|2003|) in- 
clude substantial circumstellar emission from 10 -40/im from 
old stars. This emission includes a local maximum that ap- 
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Fig. 6. — Spitzer IRAC and MIPS band flux ratios for stellar emission from various old stellar population models as a function of age (i.e. with no ISM or 
circumstellar dust contribution). The flux ratios shown are /s//3.6 (top left), /24A/3.6 (top right), fs/f 4.5 (bottom left), and /24AA.5 (bottom right). The model 
points from Maraston (2005) are for Z = 0.02 (blue circles), Z = 0.04 (blue squares), and Z = 0.07 (blue triangles). The model points from Conroy et al. (2009) 
are for Z = 0.019 (red circles) and Z = 0.035 (red squares). The empirical ratios (black solid line) and a 5000 K blackbody (dotted line) are also shown. The 
model s are a r easo nable representation of the IRAC flux ratios, but not the MIPS to IRAC flux ratios, which is evidence for a hot, circumstellar dust component. 
See §|y]and[01for details. 

4382, NGC 4570, NGC4621), and these data show clear ev- 
idence for dust emission beyound 8/im that is well fit by 
a dus ty silicate circumstellar envelope model (Br essan et al.l 
119981) . This emission spectrum is also well fit by a scaled ver- 
si on of IRS data for stars in the globular cluster 47 Tuc studied 
bv lLebzelter etaLl d2006h . 

For galaxies with a substantial amount of cold, interstellar 
dust (> 1O 6 M ), such as late-type galaxies, the circumstellar 
dust emission at 24/im is relatively negligible. However, most 
of the galaxies in our sample are substantially more dust poor 
for their stellar mass. It is consequently important to estimate 
and subtract the circumstellar contribution to the 24/im MIPS 
photometry to obtain an accurate estimate, or upper limit, for 
the cold dust component. In the next section, we describe our 
dust model fits and how we incorporate the circumstellar dust 
emission into the interstellar dust mass estimates and upper 
limits. 

5. INTERSTELLAR DUST 

Approximately half of these early-type galaxies are de- 
tected in the MIPS 70/im and 160/im bands, which indicates 
the presence of cold, interstellar dust. In this section we de- 
scribe how we fit our IRAC and MIPS photometry for these 
galaxies with the dust models of lDraine & Lil (12007b , estimate 
the total dust mass, and compare the properties of the dust in 
early-type galaxies to dust in other morphological types. We 
also derive upper limits to the cold dust mass for galaxies that 




Fig. 7. — Same as Figure [5]for the ratio of 24/im to 3.6/im (top panel) 
and to 4.5/im (bottom panel). Red points correspond to galaxies without 
evidence for dust, both based on upper limits at 160/im and the absence 
of dust structure in HST imaging. Blue points correspond to galaxies with 
160/im detections and dust structure in HST imaging. The galaxies with 
no evidence for dust have median ratios of /24/A5 = 0.252 ± 0.072 and 
fiA/h.6 =0.151 ±0.045. The errorbar in each panel represents the average 
uncertainty in the measurement of the ratio. 

proximately falls within the 24/im band and has been i den- 
tified with silicate emission in oxygen-rich AGB stars (|Suhl 
120021) . iBressan et all (120071) studied IRS Spitzer spectroscopy 
of many early-type galaxies, including a number of the in- 
active galaxies in our sample (NGC 4365, NGC 4371, NGC 
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we do not detect at 70/im and 160/im. 

All of the galaxies with dust lanes in the HST structure maps 
are also de tecte d at 8/im. We use the empirical scale factor 
derived in 34.1l to measure the PAH emission in this bandpass 
and examine the spatial distribution of the PA Hs. We also 
estim ate the fraction of the dust in PAHs with the lDraine & Lil 
( 2007) models for the galaxies with suitable data and compare 
them to other studies. 

5.1. Fit to Draine & Li Models 

Most of the dust in galaxies is in the interstellar medium. 
This dust absorbs energy from starlight, which then re-emits 
this radiation at mid- and far-infrared wavelengths. The emis- 
sion spectrum of this dust can be used to estimate the total 
mass in dust, as well as other properties, given a model for the 
typical grain size distribution and composition. These mod- 
els are used to characterize the absorption and re-emission of 
the Interstellar Radiation Field (ISRF) as a function of wave- 
length and may consequently be used to derive the dust mass 
from the observed luminosity. 

We employ the models of Draine & Li ( 2007) to character- 
ize the properties of the dust in these galaxies and derive dust 
masses. These models assume that dust is comprised of a mix- 
ture of carbonaceous and silicate grains, and varying amounts 
of PAH particles. The main observational constraints on the 
grain size distribution, relative mixture and opacities of car- 
bonaceous and silicate grains, and the properties of the PAH 
particles are the wavelength-depend ent extinction observed in 
the Milky Way (Mathi s Tet al.ll 19831), long-wavele ngth obser- 
vations of the Milky Way (Finkbeine r et al.l fl999). and stud- 
ies of the properties of PAH emission fro m the Milky Way 
and nearby galaxies (e.g. lSmith et al.ll2007l) . 

The mid- to far-infrared emission from dust grains depends 
on the incident ISRF, and thus a characterization of the dis- 
tribution of starlight intensities is also necessary to infer the 
properties of the dust. Integrated measurements of entire 
galaxies will typically include emission from dust that has 
been heated by a wide range of different values of the local 
ISRF. While some studies perform detailed radiative transfer 
calculations to estimate t he local value of the ISRF for in- 
dividual dust grains (e.g . iPopescu et al.l 120001: iGordon et al.l 
120011: iPiovan e t al. 200Q), this approach typically relies on 
many assumptions and may not be necessary to obtain rea- 
sona ble estimates of the dust mass and other dust proper- 
ties. iDraine &TH (12007b propose a simple parametrization 
of the ISRF where a fraction 7 of the dust mass is exposed 
to high-intensity radiation, such as may occur near the sites 
of star formation, and the remainder is exposed to the aver- 
age ISRF. T he dust exposed to a wide range of ISRF inten- 
sities, which Drai ne & TP (|2007) refer to as dust associated 
with Photodissociation Regions (or simply PDR dust) is ex- 
posed to a power-law distribution of intensities from U min to 
U m ax, where U is a dimensio nless scale fac t or suc h that U = 1 
for the local ISRF derived bv lMathis etaLl (119831) . U min is the 
minimum value of the ISRF, U max is the maximum value, and 
the power-law index is a. The dust exposed to the average 
ISRF is exposed to a radiation field of intensity U m m and is 
referred to as diffuse dust. The diffuse dust is expected to 
be relatively cold and dominate the emission in the 70/im 
and 160/im bands, while if present the PDR dust may be im- 
porta nt at shorter wavele ngths, most notably the 24 fim MIPS 
band. IDraine et"aT] (12007b find that the SEDs of galaxies in the 
SINGS sample are satisfactorily reproduced with U max = 10 6 
and a power-law index of a = 2. We fix these two model pa- 



rameters at those values for this study as well. The ISRF for 
the dust models is consequently just characterized by U m in and 
7. As early-type galaxies have relatively little star formation, 
we do not expect large values of 7; however, 7 may not be 
completely negligible as some early-type galaxies have mod- 
est amounts of star formation. The PDR component may also 
represent dust near the active nucleus in the small number of 
moderate luminosity AGN in this sample. 

Very small dust grains and PAHs are subject to single- 
photon heating to high temperatures and thus emit a t much 
shorter wavelengths than larger grains (Sellgren 1984). PAHs 
are also very likely to be the cause of the emission features ob- 
served at 3.3, 6.2, 7.7, 8.6, and 1 1.3 /im in the Milky Way and 
most galaxies, and in particular the 7.7/im feature appears to 
dominate the nonstellar emission in the IRAC 8/im band. This 
emission has been identified with C-C stretch an d C-H bend- 
ing m odes, primarily from ionized PAHs. The IDraine &TH 
( 2007) models include PAHs as small as 20 carbon atoms and 
note that PAHs with < 1000 carbon atoms dominate the emis- 
sion at wavelengths A < 20/xm. They present seven differ- 
ent Milky Way dust models that have different size distribu- 
tions of very small grains and different PAH fractions qpAH 
of 0.5% to 4.6%, where qpAH = 4.6% is a reasonable match 
to the Milky Way. Because each of these seven dust models 
has a different PAH fraction, q P AH uniquely describes the dust 
model. The dust emission from a given galaxy is consequently 
characterized by qpAH, U m i n , 7, and the total dust mass M Aust- 
in the next subsection we describe how we derive these four 
parameters. 

5.1.1. Flux Ratios 

The 8/im, 24/im, 70/im, and 160/im bands may be domi- 
nated by dust emission. We use the data fr om these fou r bands 
to determine the four parameters of the IDraine &TH ([2007 ) 
model: qpAH, U m i n , 7, and M dust • Our approach uses the graph- 
ical procedure described in Drai ne &Lil (1200 7) to estimate the 
best-fit dust model. This procedure employs three flux ratios 
to determine the model that is the best fit to the spectral shape. 
The total luminosity is used to estimate the dust mass in the 
context of the best-fit model. The ratios defined in that paper 
are: 
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The superscript "ns" refers to the nonstellar contribution to 
that bandpass. For the 8/im band, this is the value after sub- 
traction of the em piric al stellar contribution with the proce- 
dure described in §14.11 For the 24/im band, this is the value 
after the subtraction of the empirical stellar and circ umst ellar 
dust contribution to the 24/im band as described in 34.21 The 
measured values of these ratios are listed in Table [3] 

Not all of the galaxies with detections at 8/im, 70/im, and 
160/im have a value for Pj 9 in Table [3] because in three cases 
(NGC3377, NGC4026, and NGC5077) the flux of the non- 
stellar component at 8/im is consistent with the uncertainty 
in the empirical scale factor and our aperture measurements, 
even though weak PAH emission is apparent in Figures [2] -HI 
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Fig. 8. — Infrared spectral energy distributions and model fits. The I RAC and M IPS photometric measurements (blue circles or arrows) are shown with a 
stellar model from Maraston (2005), diffuse and PDR dust models from Draine & Li (2007) (red, dotted lines), and a simple circumstellar dust model (blue, 
dotted lines) that corresponds to a fixed fraction of the 4.5/im emission. The sum of these components is also shown (red, dashed line). The dust mass estimate, 
as well as ei ther the best-fit or assumed dust model parameters (qpAH, U m i n , and 7) are listed below each galaxy's name. The dust model fits are described in 
section 35.11 



For a similar reason, no P24 value is quoted for seven galax- 
ies (the previous three cases, as well as NGC2768, NGC2787, 
NGC4278, and NGC4589) because there is more substantial 
variation in the empirical scale factor that corrects for the con- 
tribution from (mostly) hot, circumstellar dust in this pass- 
band. 

The ratio P7.9 is effectively a measure of the fraction of the 
dust emission that is radiated by PAHs, and thus provides the 
strongest constraint on the dust model and the qpAH parameter. 



iDraine &~Q (|2007) note that because the 8/im emission is al- 
most exclusively due to single photon heating, P7.9 is not very 
sensitive to the starlight intensity. The exception is that if the 
PAH fraction is very small, exposure to a more intense ISRF 
(larger 7) will lead to enhanced PAH emission that may pro- 
duce values of Pj g that are comparable to models with larger 
qPAH and smaller values of 7. P7.9 is less sensitive to 7 for 
larger values of qpAH because these galaxies exhibit substan- 
tial PAH emission even when 7 = 0. 
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Fig. 9. — Continuation of Figure [8] 
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Fig. 10. — Continuation of Figure[8] 
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Fig. 11. — Continuation of Figured] 



The value of 7 is mostly determined by the ratio P24. This 
ratio is sensitive to the fraction of the dust emission near 
24/im and is dominated by hot, P DR dust (7 > 0) if a h ot com- 
ponent is present. According to iDraine &TH (120071) . when 
larger grains are heated by starlight intensities above U ~ 20, 
they will contribute to the 24/im emission. In contrast, dust 
that is heated by starlight intensities ofU ~ 0.1 - 10 is suffi- 
ciently cold that it emits most of its luminosity at longer wave- 
lengths, such as the 70/mi and 160/im bands. 

The coldest dust component is diffusely distributed in the 
ISM and heated by the average ISRF of the galaxy, U m i n . The 
ratio R71 is sensitive to U m i n because the intensity of the ISRF 
is correlated with the characteristic temperature of the cold 
component, and this affects the ratio of the flux in the 70/iin 
and 160/im bands. For small values of 7 (on order a few per- 
cent or less), the ratio Rj\ is a very sensitive to JJ m [ n because 
the 70/im flux has a negligible contribution from PDR dust. 
For larger values of 7, the PDR dust contribution to the 70/im 
emission is more important. In galaxies with substantial star 
formation, the values of both Rq\ and P24 are larger and Rq\ is 
less sensitive to U m i n . 



TABLE 3 
Dust Model Fits 



Name 


P7.9 


P24 


R71 


qpAH [%] 


Umin 


7 


log M dust 


NGC0315 


0.13 


0.48 


1.47 


2.50 


2.0 


0.2 


6.7 


NGC0821 








2.50 


7.0 


0.0 


<4.7 


NGC1023 








2.50 


7.0 


0.0 


<3.9 


NGC2300 








2.50 


7.0 


0.0 


<4.9 


NGC2768 


0.06 




2.31 


1.12 


15.0 


0.0 


5.5 


NGC2787 


0.11 




2.52 


1.77 


15.0 


0.0 


5.2 


NGC3226 








2.50 


7.0 


0.0 


5.5 


NGC3377 






0.75 


2.50 


2.0 


0.0 


4.7 


NGC3414 








2.50 


5.0 


0.0 


5.6 


NGC3607 


0.18 


0.07 


1.53 


3.19 


7.0 


0.0 


6.3 


NGC3640 








2.50 


7.0 


0.0 


<4.6 


NGC3945 








2.50 


7.0 


0.0 


5.4 


NGC3998 


0.25 


0.44 


1.85 


3.90 


7.0 


0.15 


5.7 


NGC4026 






1.18 


2.50 


5.0 


0.0 


5.3 


NGC4138 


0.19 


0.13 


1.30 


3.19 


5.0 


0.02 


6.3 


NGC4278 


0.10 




2.09 


1.77 


5.0 


0.0 


5.4 


NGC4291 








2.50 


7.0 


0.0 


<5.0 


NGC4293 


aio 


a 19 


1.65 


1.77 


7.0 


0.04 


6.5 


NGC4365 








2.50 


7.0 


0.0 


<4.2 


NGC4371 








2.50 


7.0 


0.0 


<4.3 


NGC4382 








2.50 


7.0 


0.0 


<4.7 


NGC4406 








2.50 


7.0 


0.0 


<4.5 


NGC4526 


0J2 


o!o7 


235 


2.50 


15.0 


0.0 


6.4 


NGC4550 








2.50 


7.0 


0.0 


4.9 


NGC4570 








2.50 


7.0 


0.0 


<4.2 


NGC4578 








2.50 


7.0 


0.0 


<4.3 


NGC4589 


0.08 




1.37 


1.12 


3.0 


0.0 


6.1 


NGC4612 








2.50 


7.0 


0.0 


<4.2 


NGC4621 








2.50 


7.0 


0.0 


<4.1 


NGC4636 








2.50 


7.0 


0.0 


4.7 


NGC4649 








2.50 


7.0 


0.0 


<4.2 


NGC4694 


0.23 


0.14 


1.40 


3.90 


5.0 


0.02 


6.1 


NGC5077 






1.95 


2.50 


12.0 


0.0 


5.7 


NGC5273 


0.15 


(X28 


2.45 


3.19 


15.0 


0.06 


5.4 


NGC5308 








2.50 


7.0 


0.0 


<4.6 


NGC5557 








2.50 


7.0 


0.0 


<4.9 


NGC5576 








2.50 


7.0 


0.0 


<4.6 


NGC7619 








2.50 


7.0 


0.0 


<5.1 



Note. — Dust model results of the galaxies in the sample. Columns are: (1) 
Galaxy name; (2) P7.9; (3) P24; (4) R71; (5) q PAH \ (6) U min ; (7) 7; (8) log of M dmt 
in solar masses. Galaxies with measurements of P7.9, P24 and/or R71 were fit with 
the Draine & Li 1 2007) models. The remaining galaxies were fit with a fixed set of 
parameters: (qpAH, Umin, 7) = (2.50%, 7.0, 0.0). Dust mass upper limits are for this 
fixed model and the 3a upper limits at 70/im and 160/im. More details are provided 
in §5.1. 



5.1.2. Dust Model Parameters 

There are eight gala xies in our sample with sufficient data 
to constrain all three IDraine &TH (|2007) model parameters. 
For these g alaxies we fo l low th e graphical procedure recom- 
mended by IDraine &TH (12007b and first estimate qpAH from 
the values of P7.9 and Rj\, and then use the values of P24 and 
Rji to estimate 7 and U m [ n . Four additional galaxies have just 
P7.9 and 7?7i, which is sufficient to constrain the dust model. 
A histogram of qpAH for these 1 1 galaxies is shown in the first 
panel of Figure [T2l The median value of qpAH is 2.50% for 
these galaxies. The second panel shows the distribution of 7 
for the subset of eight galaxies with sufficient coverage of the 
diffuse emission. The median value value of 7 for these eight 
is 0.03. Finally, there are three galaxies with just Rq\. For 
these galaxies we assume qpAH = 2.50% and 7 = (to provide 
a more conservative constraint on the dust mass than 7 = 0.03) 
and then estimate U m i n . The median value of JJ m { n is 7.0. 

With the dust model fixed, the dust mass is: 

M dust = -^ ((vfv)24 + (vfv)n + (vfv)i6o)D 2 . (4) 

The quantity ^ is uniquely determined by the three parame- 
ters qpAH, 7, and U m i n that describe the dust model, the quan- 
tity (U) accounts for starlight heating in both PDRs and the 
diffuse ISM, and D is the distance. Dust mass histograms are 
shown for different morphological types in Figure [13] and for 
active and inactive galaxies in Figure fT4l 

The distributions of qpAH, 7, and U m j n are generally differ- 
ent from the distributions found by Dra ine et al.l (2007) for 
the SINGS sample. The qpAH values are generally lower, 
with the median of Qpah = 2- 5 % bel ow the median value of 
3.4% found by IDraine et all (120071) for the 61 galaxies in 
the SINGS sample with similar IRAC and MIP S data. The 
8 PAH feature is mostly due to ionized PAHs (Jobli neFaT] 
1994), so the weakness of this feature may indicate a greater 
fraction of neutral P AHs, rather than a smaller PAH fraction. 
iKaneda et al.1 d2Q05|) proposed this explanation based on their 
study of IRS spectra for four elliptical galaxies, which showed 
the short- wavelength PAH features were weak or absent, but 
the longer wavelength features at 11.3/zm and 12. 7 /mi were 
prominent. Their sample includes NGC4589, for which our 
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Fig. 12. — Distribution of the dust model parameters qpAH (left), 7 (middle), and U m in (right) in galaxies with sufficient detections to constrain one or more of 
these quantities. Note that NGC3998 has 7 = 1 and consequently does not fall within the range shown in the middle panel. 




Fig. 13. — Distribution of dust masses and upper limits as a function of 
morphological type. 
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Fig. 14. — Distribution of dust masses and upper limits as a function of 
nuclear activity. 

best model fit yields qpAH = 1.12%, yet the longer wave- 
length PAH f eatures are clearly present in the IRS spectrum. 
iKaneda et al.1 d2005l) suggest the relative absence of UV radia- 
tion in ellipticals could produce physical conditions amenable 
to primarily neutral PAHs. Such conditions are also indicated 
by the spectral ene rgy distribution of the bulge of M31 by 
iGroveset all (120121) . 

Other potential explanations include the preferential de- 
struction of small particles, such as due to sputtering in hot 
gas or an AGN, and low gas-phase metallicities. The domi- 
nance of larger PAHs ov er smaller ones, which could be a con- 
sequence of sputtering ( Schutte et al. 1993; Tsai & Mathews 



119951: iMicelotta etail 120101) . would weaken the 8/im fea- 
ture. The qpAH value is also observed t o decline in galaxies 
with low gas-phase oxygen abundance. iDraine et"aTT (2007) 
find that the median PAH fraction is 1.0% for galaxies with 
gas-phase oxygen abundances below 12 + log 10 O/H < 8.1 (it 
is also low for a small fraction of higher metallicity galax- 
ies). Even though our sample is expected to possess solar to 
super- solar stellar abundances, the cold gas abundance may 
be quite low if it is the remant of an accreted, low-mass satel- 
lite. Finally, relatively luminous AG N have also be en impli- 
cated in the destruction of PAHs (iLutz et al.l [20081) and all 
but three of the galaxies with detections in all seven bands 
are AG N, although they are low luminosity AGN. Smit h et aTl 
(2007) found that AGN may modify the relative distribution 
of intensities of various PAH features, and in particular sup- 
press those from 5 - 8/im. However, low-luminosity AGN 
do not appear to n egatively effect the PAH abundance, as 
Draine et al. (2007) found that the median q P AH fraction was 
actually slightly higher (3.8%) for AGN than the median of 
their full sample, as well the median of 28 galaxies with HII 
nuclei (3.1%). 

Nearly all of the galaxies have 7 < 0.1 and the median 
value is 7 = 0.03. This is consistent with the near- absence 
of star formation in these galaxies. The small values of 7 
also mean that the model fits and dust masses are more re- 
liable, as in this regime the parameter Rji provides a bet- 
ter estimate of U m { n . The median U m i n value for this sample 
is 7.0, which is several times higher than the SINGS sam- 
ple median of 1.5. The median for the subsample of SINGS 
galaxies with SCUBA data is in better agreement with our 
data. In the case of SINGS galaxies with SCUBA obser- 
vations, the cold and PDR dust components are better sep- 
arated, and thus U m i n is better measured for those galaxies. 
We may be able to measure larger values of U m in without 
longer-wavelength sub-mm observations because our galax- 
ies have lower values of 7. Early-type galaxies may also 
have more intense ISRFs than spiral galaxies due to the pres- 
ence of hot, evolved stars and higher stellar densities. This 
is also su ggested by th eir hi gher dust temperat ures relative 
to spirals dBendo > et al.M2Q13b ISkibba et alJl20lU ISmith et al.1 
2012: lAuld et al] 12013). Emission from hot, evolved stars 
contributes to and may dominate the L INER emission spectra 
from some galaxies dSarzi et alll2010|) . Longer wavelength 
observations would help further constrain the value of U m m, 
while far-infrared line ratios would hel p to constrain the spec - 
tral shape of the ionizing photons fe.g. lMalhotra et al. 2000). 

An important caveat that affects the interpretation of the 
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dust model parameters is that uncertainties, or intrinsic vari- 
ations, in the empirical scale factors derived in Section 14.11 
increase the uncertainties in the dust model parameters. As 
noted previously, the data shown in Figures [5] and [7] exhibits 
evidence for intrinsic variation in these quantities, and some 
variation is expected due to stellar population differences. 
This is particularly relevant for P7.9 and which largely 
determine qpAH and 7, and so these parameters are not as well 
constrained as U m i n . In practice, as 7 is expected to be small 
due to the general absence of star formation in these galaxies, 
the biggest potential uncertainty lies with qpAH- 

The measured dust masses for our sample are Mdust ~ 
10 5-6 5 Mq, which are two or three orders of magnitude lower 
than the dust masses for spiral galaxi es of comparable size, 
but similar to other early-type galaxie s (Goudfr ooij & de Jongl 
119951: iDraine et al.ll2007b ISmith et al.l[2012l) . The dust masses 
that result from this procedure are fairly robust, even though 
they are only based on measurements out to 160/im. This 
is true even if the model parameters qpAH, 7, and U m i n are 
not well constrained by the data because qpAH has a negli- 
gible impact on the dust mass estimate and 7 = is a good 
approximation f or most galax i es. Fo r the full range of mod- 
els presented in IDraine &TH (|2007), we also note that the 
value of ^ that is set by the best-fit model only varies from 
~ 0.044-0.066 g (erg s -1 ) -1 . This indicates that even if 
the dust model were not constrained at all by the data and 
one simply adopted the average of this range, there would 
only be a ^ 20% uncertainty in the dust mass within the 
context of these models. An additional uncertainty is the 
ISRF because the absence of flux measurements at wave- 
lengths longer than 160/im do not constrain a potentially sub- 
stantial mass of cold dust. However, the available data can 
constrain the ISRF quite well because these galaxies have 
little star formation. Another estimate of the uncertainties 
comes from the overlap between several g alaxies in our study 
and previ ous investigatio ns with ISO (Xilou riset al.l 12004). 
Spitzer (Temietal. 200 3), and longer- w avelength Herschel 
data (ISmith et al.ll2Q12h lAuld et al.1 (20131) . The dust masses 
derived by those previous studies assumed a modified black 
body and mostly agree with our estimates to within a factor 
of two. The notew orthy exception is NGC3945, for which 
Smit h et aTl (|2012) derive over an order of magnitude larger 
dust mass. The difference could be ascribed to the fact that 
we use a smaller aperture, although our 70/im flux is within 
5% of the IRAS 60/im value. However, due to the large uncer- 
tainty in this IRAS measurement it is still possible that we are 
missing extended emission from dust in this object. We may 
be able to put tighter constraints on the dust mass for this ob- 
ject with future analysis of Herschel photometry. We also note 
we were unable to obtain a good measurement of NGC3945 at 
160/im. We do retain NGC3 945 in our sample. Another com- 
parison was performed by Gor don et al.l (120101) . who found 
that dust masses estimated from < 160/im data can differ by 
10-36% relative to dust mass estimates that include Herschel 
measurements at 250/im and 350/im, although this study was 
based only on the observations of the Large Magellanic Cloud 
and therefore does not include uncertainties due to differences 
in grain size and ISRF intensity distributions . Nev ertheless, 
this study and a recent one by lAniano et al.l (1201 2|) of two, 
nearby spirals that also compared Spitzer and Herschel pro- 
vides good evidence against substantial masses of cold dust. 

A final uncertainty is that the physical dust model that 
works so well for the Milky Way may not be appropriate 



for early-type galaxies, that is the size and composition are 
different because of different formation and processing histo- 
ries. For example, sputtering can reduce the fraction of small 
grains relative to the initial grain size distribution. The gen- 
erally smaller values of qpAH in early-type galaxies relative to 
spirals may indicate that this is the case, although there is also 
allowance for this difference in the Draine & Li (2007) mod- 
els. The study of four ellipticals by Kaneda et al. (2005) also 
provides some evidence that the shape of the radiation field is 
softer in ellipticals than in the Milky Way, which may not be 
adequately compensated by an adjustment in U m i n . Provided 
all early-type galaxies are different in the same manner, then 
the relative dust masses we have derived here should never- 
theless be robust. The true uncertainties in the dust masses 
are likely dominated by these model details and are difficult 
to quantify. Under the assumption that the dust model is a rea- 
sonably good match to early-type galaxies, we estimate that 
the uncertainties in the dust masses are on order a factor of 
two or three due to a combination of aperture differences, cal- 
ibration uncertainties, and differences (or uncertainties) in the 
dust models. 

5.1.3. Upper Limits 

Approximately half of the galaxies in our sample are not 
detected at 70/im and 160/im. These non-detections never- 
theless provide very useful constraints on the total dust mass 
in these galaxies because of the extraordinary sensitivity of 
the Spitzer MIPS instrument. To calculate upper limits on the 
dust masses in these galaxies, we employ a dust model with 
qPAH = 2.50%, 7 = 0, and U m i n = 7.0. The values of qpAH and 
U m i n are representative of the galaxies with the best data (see 
Figure [12]). We chose to use 7 = 0, rather than the median 
value of 7 = 0.03, because this leads to a more conservative 
constraint on the dust mass. 

We use this dust model to calculate the largest value of 
the dust mass that is consistent with the 3a upper limits at 
70/im and 160/im. These upper limits are listed in Table [3] 
and shown in Figures Q3] and [14] The figures clearly show 
that the upper limits on the dust masses are approximately an 
order of magnitude lower than the measured dust masses. 

All of these upper limits are for galaxies that show no ev- 
idence of dust lanes within 10 0s of parsecs of their n uclei 
based on our earlier HST study ( Simoes Lopes et al. 20071 and 
also shown in the leftmost panels of Figures [2]- [4]). While the 
structure map technique is insensitive to uniformly distributed 
or very diffuse dust, the Spitzer data place very strong upper 
limits on the mass of such dust. The PAH images are also 
valuable in this respect, as we describe in the next subsection. 

5.2. PAHs 

As a second search for diffuse, extended dust that may not 
have been present in our HST study, we created PAH images 
for all of the galaxies. These images are shown in the mid- 
dle panels of Figure [2] - [4] The PAH images are the 8/im 
image minus a scaled version of the 4.5/im image, where the 
scale factor accounts for the relative stellar emission in the 
tw o ima ges. The calculation of the scale factor is described 
in 34.11 The pair of images were convolved to the same PSF 
prior to subtraction. All of the galaxies with dust lanes in the 
HST structure maps also have clear PAH emission, with the 
exception of NGC4636, which has only marginally-detected 
dust lanes in the HST data. Conversely, none of the galaxies 
with no evidence for dust at visible wavelengths exhibit clear 
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evidence of PAH emission. A number of the PAH images of 
the galaxies exhibit artifacts at their centers. These are consis- 
tent with artifacts of the stellar continuum subtraction, which 
could be caused by uncertainties in the scale factor, the PSF 
shape, and centroid errors. 

There is a good correspondence between the morphology 
of the dust lanes in the HST data and the distribution of PAH 
emission for the small subset of galaxies where the HST data 
are of sufficient quality and the dust lanes sufficiently ex- 
tended. Some particularly striking examples are NGC4138 
and NGC4293. In a few cases the structure maps show sub- 
stantial asymmetries in the dust distribution that are not mir- 
rored in the PAH images. The most likely explanation is 
that the dust lanes viewed in absorption are substantially less 
prominent on the si de of the gala xy that is tilted away from 
our perspective (e.g. Hubble 19431). 

5.3. Distribution of Dust Masses 

The dust mass distribution shown in Figure [13] shows that 
most detections have 10 5-6 M of dust. There are only 
three detections below 10 5 M (NGC3377, NGC4550, and 
NGC4626), while there are many upper limits that extend an 
order of magnitude lower and suggest a relative absence of 
galaxies with < 10 5 M of dust. For example, if there were 
as many galaxies with dust masses in the range 10 4-5 M as 
in the range 10 5-6 , we would expect ten galaxies in this lower 
range and that we would detect half of them, yet only one is 
detected. This apparent dichotomy is also supported by the 
similar dichotomy in the detection of dust lanes in HST im- 
ages of these galaxies (see Table 1), as we do not expect the 
Spitzer and HST images to have comparable dust mass sensi- 
tivity. 

One possible explanation is a difference in morphological 
type. Figure Q3] splits the detections and upper limits by mor- 
phology into ellipticals, SOs, and E/SOs. We see no evidence 
that true ellipticals or true SO galaxies are more or less likely 
to have interstellar dust. The sensitivity of the datasets for 
each morphological type are also comparable. The upper lim- 
its on the dust mass in undetected ellipticals and lenticulars 
both span the full range of upper limits from approximately 
10 4-5 M . We used th e Astronomical Survival A nalysis 
(ASURV) programs from Feige lson & Nelson! (|1985|) . as im- 
plemented in the I RAF STSDAS STATISTICS package, 
to compare the elliptical and lenticular subsamples. We ran 
five different tests through the twosampt task and these in- 
dicated that the distributions are consisten t. These simi l ar dis- 
tributions are in contrast to the results of Smi th et ail (j2012) 
with the Herschel Reference Survey, who detected dust in 
24% of ellipticals and 62% of SOs with SPIRE. This differ- 
ence in detection rate between the two morphological types 
could be due to their larger sample o f 62 early-type gal axies. 
Another potential explanation is that lCortese et al.l (|2012|) and 
found that the dust to stellar mass ratio is lower for ellipticals 
than lenticulars. If our ellipticals were more massive than our 
lenticulars, this would also explain their similar dust masses. 

We investigated if there is a correlation between dust de- 
tections or upper limits and cluster membership. Environ- 
ment may correlate with the absence of dust because the hot 
atmospheres of clusters may lead to an increase in the rate 
of dust destruction, although many cluster early-type galax - 
ies do maintain their own atmospheres (e.g. ISun et al.l l2005). 
With regard to environment, clusters of galaxies have substan- 
tially fewer gas-rich dwarf galaxies and sufficiently high rel- 




FiG. 15. — Dust mass M^ust vs - luminosity in at 4.5/mi L4.5 for galaxies 
with good dust mass estimates (solid points) and upper limits (arrows). Ac- 
tive galaxies are marked with red symbols and inactive galaxies with blue 
symbols. The upper limits are based on the fiducial dust model described in 
§5.1.3 and the 3a upper limits on the flux at 70/im and 160/xm. 

ative velocities that mergers are more rare. We searched for 
differences as a function of environment, and whether or not 
environment could explain the dust ditchotomy, with the 14 
Virgo members and 27 nonmembers in our sample and de- 
tect dust in 18/27 (66.6%) nonmembers, compared to 5/14 
(45.4%) Virgo members. This difference suggests a slightly 
larger incidence of dust in galaxi es outside of Virg o, although 
it is not statistically significant. Smit h et aT] (120121) performed 
a similar analysis for the early-type galaxies in the Herschel 
Reference Survey and found a similar trend, but that it was 
also not statistically significant. We also investigated whether 
or not membership influenced the detection of dust in the in- 
active sample. We detect dust in 2 of the 1 1 inactive Virgo 
members, compared to three out of 1 1 of the inactive galaxies 
outside of Virgo. There is thus no evidence that membership 
in the Virgo cluster is correlated with a lower incidence of 
dust in inactive galaxies. 

We also compared the dust detection rate with stellar lumi- 
nosity, which correlates well with stellar mass. More massive 
galaxies generally have higher ISM temperatures, which will 
lead to faster dust destruction rates. This naively suggests an 
anticorrelation between dust mass and stellar luminosity, al- 
though no such anticorrelation is present in our data. In Fig- 
ure[T5]we show the dust masses and upper limits as a function 
of L4.5, a proxy for stellar mass. While the minimum detected 
dust mass is a function of stellar luminosity, this is because 
more luminous galaxies are more rare and thus at larger dis- 
tances. 

The dichotomy between the dust mass estimates and up- 
per limits thus suggests that there are distinct populations of 
dusty and non-dusty early-type galaxies, rather than a con- 
tinuum of dust masses. To explore this possibility further, we 
resampled our data many times in order to determine if the ap- 
parent dichotomy is due to the sensitivity distribution of our 
data, rather than an intrinsically nonuniform distribution of 
dust masses. To perform these calculations, we estimated the 
minimum flux sensitivity of the 160/im data for each galaxy. 
We then generated a uniform distribution in log Md ust , ran- 
domly assigned dust masses from this distribution to galaxies 
in our sample, calculated the flux at 160/im based on a ran- 
domly selected galaxy's distance for the dust model used to 
estimate upper limits, and determined if this flux would corre- 
spond to a detection or upper limit based on the sensitivity of 
the same galaxy's dataset. We generated 1000 realizations of 
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our dataset and found that the sensitivity and distance distribu- 
tion of our dataset does not produce the observed dichotomy 
between detections and upper limits for an intrinsically uni- 
form distribution in log Mdust approximately ~ 95% of the 
time. 

While this is formally marginally significant, an important 
caveat to this statistical analysis is that it does not include vari- 
ations due to the best-fit dust model, which may add a factor 
of two uncertainty to the dust mass, the possibility that we 
may systematically underestimate the dust mass by tens of 
percent in some cases, particularly for nearby galaxies with 
extended emission, and potential correlations with morphol- 
ogy, stellar mass, and environment. We plan to reexamine this 
potential dichotomy with a larger sample. 

5.4. Correlation with AGN 

The distribution of dust masses is strikingly different for ac- 
tive and inactive galaxies. Galaxies classi fied as ac t ive, tha t 
is either LINERs or Seyferts according to iHo et al.l (Il997al) . 
are all detected and have dust masses that range from M = 
10 5-6 5 M . In contrast, only four galaxies classified as in- 
active are detected and the remaining 16 have upper lim- 
its that range from M ~ 10 4-5 M©. This result is similar 
to the strong correlation between activity and the presence 
of dus t structure within the cen tral kpc in our earlier HST 
study (ISimoes Lopes et al"1 l2007). yet these new observations 
demonstrate that the galaxies without dust lanes have at least 
one to two orders of magnitude less dust than those with de- 
tections. As in the previous section, we used the ASURV pro- 
grams from iFeigelson & Nelsonl (11985b to compare the dust 
mass distributions of the active and inactive samples. All 
five tests available through the twosampt task returned zero 
probability that the dust masses were drawn from the same 
distribution. 

There has been substantial discussion in the literature of the 
extent to which the emission lines characteristic of LINERs 
are actually powered by accretion onto supermassive black 
holes. Other proposed candidates are fast shocks, interac- 
tion with the hot phase of the ISM, and photoionization by 
UV emissio n from old s tellar populations (see lHoll2008l for 
a review). ISarzi et ail (|2010|) performed a detailed analy- 
sis of emission-line ratios for many early-type galaxies and 
concluded that the emission-line luminosities of only a mi- 
nority of LINERs are dominated by accretion. Nevertheless, 
AGN have been confirmed in the majority of LINERs based 
on measurements of radio cores, X-ray cores, and the sur- 
face brightness profiles of various emis s ion lines relative t o 
the stellar continuum iFilho et al I (|2004|); iNagar et al.l d2005|); 
Filho et al.1 (120061) ; iFlohic et al.l (12006b : lHd " d20Q8h : ISarzi etall 
(2010). Some fraction of the cold ISM in the majority of these 
dusty early-type galaxies is therefore fueling their central, su- 
permassive black holes, even if their emission line luminosi- 
ties are not dominated by nuclear accretion. We hypothesize 
that the correspondence between active galaxies and dust in 
our sample is simply due to the presence of a sufficient quan- 
tity of cold material, some of which has accreted onto the 
black hole, and some of which is ionized and produces the 
visible-wavelength emission lines, rather than a causal rela- 
tionship between the presence of dust and AGN. 

6. ORIGIN OF THE DUST 

The main motivation for this study was to use dust mass es- 
timates for early-type galaxies to constrain the origin of their 



dust. Our estimates demonstrate that dusty early-type galax- 
ies have 10 5-6 - 5 M of dust, while the remainder have at least 
one order of magnitude less. Our sample is representative of 
the early-type galaxy population as a whole because the orig- 
inal selection of th e sample w as from the magnitude-limited 
Palomar Survey of IHo et all (Il995|) . In ISimoes Lopes etall 
d2007h . we selected an equal number of active and inactive 
galaxies from all of the early-type galaxies in the Palomar 
survey that had visible- wavelength HST images. Half of this 
sample were active galaxies and half were not, which is equiv- 
alent to the relative frequency of active and inactive galax- 
ies in the early-type galaxy population. Sim oes Lopes et al.l 
(2007) then found that all of the galaxies classified as active 
had dust lanes within the central kpc, whereas only 25% of 
the inactive galaxies had them. While we have more inac- 
tive galaxies than active galaxies in this sample, we still find 
that all active galaxies have interstellar dust and find a sim- 
ilarly small fraction (5/22 or 23%) of the inactive galaxies 
have dust. We therefore conclude that ~ 60% of all early-type 
galaxies have dust, including all of those classified as active, 
and ~ 40% do not. 

The presence of dust in any early-type galaxy has long 
been a puzzle because thermal sputtering of dust grains in 
the hot (T > 10 6 K) gas that comprises the bulk of their 
ISM should destroy dust o n the timescale Tdust ~ 2 x 10 4 yr 
(Draine & Salpeter 1979b). These galaxies clearly do have 
evolved stars that produce dust, as indicated by the strong evi- 
dence for hot, circumstellar dust, and this dust must be ejected 
into the ISM at an approximately constant rate due to the age 
of their stellar populations. The expected rate of mass loss, 
the dust-to-gas ratio of this material, and the dust destruction 
timescale, can be used to estimate the steady-state dust mass 
of these galaxies. As noted in the Introduction, for reasonable 
estimates of the mass-loss rate, dust-to-gas ratio, and dust de- 
struction time in hot gas, the steady-state dust mass should 
be 10- 100(r^/2 x 10 4 yr) M , that is none of the galaxies 
should have been detected. 

The dust destruction timescale could be longer if the tem- 
perature of the hot phase is closer to 10 5 K, or if some of 
the gas and dust cools rapidl y and settles into a cooler phase 
(Mat hews & Brighentil 120031) . although the dust destruction 
timescale needs to be many orders of magnitude longer to pro- 
duce dust masses of 10 5 M and above. A further complica- 
tion is that this calculation predicts the expected steady-state 
dust mass for all early-type galaxies. As early-type galaxies 
form a fairly homogeneous population, all early-type galax- 
ies should have similar dust production rates and destruction 
timescales. The internal origin model is consequently incon- 
sistent with the order of magnitude differences in dust mass 
for similar galaxies, such as is illustrated by Figure IT51 

The many galaxies with upper limits on the dust mass pro- 
vide an upper limit on the typical dust destruction timescale 
in the hot ISM through an inversion of the argument presented 
above. For typical early-type galaxies, the most stringent up- 
per limits on the dust mass are approximately Mdust < 10 4-4 5 . 
For a continual dust production rate of 0.1- 1 x 0.005 = 
5-50 x 10~ 4 M yr -1 , the upper limit on the dust destruc- 
tion timescale is: 

~ in7 ( Mdust,lim \ ( 0.1 \ /0.005\ 

Ti " <2x,o [—)[i^j{—) yr <5> 

where Mdust.Um is the upper limit on the dust mass, M gas is 
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the gas mass loss rate of the evolved stellar population, and 
€d g is the dust to gas ratio. This upper limit is comparable 
to the value of Tdust < 46 ± 25 Myr derived by Clemens et al. 
(1201 Q|) with a s i milar argument, as well as more recent work 
by ISmith et all (|2012|) . Both values are consistent with, al- 
though se veral orders of magnitude greater than, the calcu- 
lations of iDraine & Salp eter (1979b), and inconsistent with 
an internal origin for the one to two orders of magnitude 
more dust present in approximately half of the sample. As 
noted in the introduction, a number of previous studies have 
also pointed out the difficulties with a purely in ternal origin 
(Sim oes Lopes et afll200 7: Rowlands et al.1 l2012|) . In particu- 
lar, the recent study by Smi th et al.l (1201 2|) reached a similar 
conclusion with a similar analysis of 62 early-type galaxies 
observed by Herschel. 

The alternate, and favored, hypothesis for the origin of 
the dust is the accretion of gas-rich satelli tes (e.g. lTran et al.l 
120011: IVerdoes Kleijn & de Zeeuwl [2005). The morphol- 
ogy of the dust lanes observed in HST images indicate 
that many of them are chaotic and clumpy amd this ap- 
pears inco nsistent with steadily accumulation from stellar 
mass loss (Ivan Dokkum & Franxl 119951: ISimoes Lopes et all 
120071) . Furthermore, when smooth dust disks are present, 
they do not always share the maj or axis of the host galaxy . 
Kinematic observations of atomic (IBertola et al1 [T984. 1992; 
Sarzi et al. 2006; Morsanti et al. 2006, e.g.) and molecu- 
lar dCombes et al. 2007; Davis e t al.ll201l)) gas in early-type 
galaxies also indicate that the gas kinematics are often not 
aligned with the kinematics of the stars, which further sug- 
gests an external origin. 

The key challenges to the external origin hypothesis are that 
the merger rate of gas-rich satellites is low and the dust de- 
struction timescale is short, yet dust is present in slightly more 
than half of all early-type galaxies. That is, the external origin 
model must satisfy the constraint: 

fdust — ^-mergTdust (6) 

where fdust = 0.6 is the fraction of the population with dust 
above some minimum mass, 3l mer g is the merger rate of gas- 
rich satellites that could supply sufficient dust, and Tdust is 
the dust destruction time. Dwarf galaxies with 10 5-6 M 
of dust (comparable in size to the Magellanic clouds) would 
correspond to 1 : 300- 1 : 100 mass-ratio mergers for the 
galaxies in this sam ple. Based on the equations provided by 
StewartetalJ(2009), we calculate that the cumulative merger 
rate for mass ratios from equal mass mergers to 1 : 300 merg- 
ers is ^Stmerg = 0.07-0.2 Gyr -1 , which is more than two or- 
ders of magnitude higher than our upper limit of Tdust < 0.02 
Gyr and more than f our orders of magnitude g reater than the 
theoretical estimate (IDraine & Salpeter 1979b). These values 
predict f dust < 0.0014-0.004, or a discrepancy of 150-430. 

To evaluate the significance of this discrepancy, we esti- 
mate the uncertainty in both quantites. First, the merger rate 
we have adopted may be an overestimate because the theo- 
retical calculations predict the merger rates for all galaxies 
as a function of stellar mass, yet not all satellites galaxies 
may have sufficient cold gas. It is also possible that the the- 
oretical estimates underestimate the true merger rates. Ob- 
servationa l evide nce for an underestimate was presented by 
lLotz et all (|2011|) . who found that the observed merger rates 
ar e approximately a fa ctor of five higher than the predictions 
of IStewart et al.l (|2009|) and other recent, theoretical analyses. 
lLotz et all (1201 ll) posit that the disagreement between theory 



and observation could be because the visibility timescale for 
the mergers has been overestimated, or the theoretical rate has 
been underestimated. While the merger visibility timescale is 
arguably more uncertain that the merger rates from simula- 
tions, this comparison indicates the merger rates are unlikely 
more than a factor of a few larger than the theoretical esti- 
mates. There is consequently insufficient uncertainty in the 
merger rates to resolve this discrepancy. 

The other potential resolution is that the dust destruction 
timescale is actuall y substantially higher than both the the- 
oretical estimate of IDraine & Salpeter (1979b) and our em- 
pirical upper limit. While this requires an increase of about 
two orders of magnitude relative to the upper limit, the 
dust destruction timescale is much more uncertain than the 
merger rate. One reference value is the dust destruction 
timescale for the Milky Way Tdust ,mw = 0.4 Gyr, which is dom- 
inated by the impact of supernova shocks on the cold, neu- 
tral medium dBarlow & Silldll977b IDraine & Salpete3ll979al: 
iDwek & S calo 1980l: lJones et al.ll994l) . As the supernovarate 
in early-t ype galaxies is several times lower than in the Milky 
Way (e.g. iLi et al.ll2011|) . the dust destruction timescale may 
be several times larger in their cold ISM and plausibly on 
the order of a Gyr, provided the accreted dust does not mix 
with the hot ISM. This may be the case if a much larger dust 
destruction timescale, combined with the highest estimated 
merger rates, predicts a dusty fraction that is only a factor of 
a few below what is needed to reproduce the observations. 
This is plausibly within the uncertainties of these order of 
magnitude estimates, particularly if dust destruction is very 
inefficient in the cold ISM of early-type galaxies and the cold 
ISM does not mix efficiently with the hot phase. However, 
it is also plausible that the lifetime will be less because there 
is less cold ISM. A more precise estimate could be obtained 
with an analysis of the fraction of the supernova energy that 
impacts the cold ISM. 

While it may be possible for purely external accretion 
to work, we instead propose that the solution is a hybrid 
of internal production and external accretion. When early- 
type galaxies accrete a gas-rich satellite, they accrete sub- 
stantial amounts of both cold gas and dust. While dust 
will be destroyed by a combination of supernova shocks 
and sputtering, this destruction could be balanced by con- 
tinued dust grain growth in the accreted cold ISM. While 
this process is not constrained by observations of early- 
type galaxies, many studies of the Milky Way have con- 
cluded that > 90% of dust formation by mass occurs in 
the cold ISM j Draine & Sal pete3fl979al:lDwek & Scaloll 19801: 
lMcKedfl98^ : Draine 2009), rather than from stellar sources 
such as mass loss and supernova alone (|Barlow et al.l 1201 (A 
iMatsuura et~aT1l201 lblPunne et al.ll201 lb . The Milky Way has 
M gas = 5 x 10 9 M of cold gas that produces M dus t = 2.5 x 10 7 
M per dust destruction timescale, which corresponds to a 
steady-state dust mass of €d g M gas . While the dust destruc- 
tion timescale is likely different from the Milky Way value, 
the steady-state dust mass is likely to be within an order of 
magnitude of that predicted by the product of the Milky Way 
dust-to-gas ratio and the observed cold gas supply. 

Continued dust production in the accreted cold gas removes 
the main weakness of the internal origin hypothesis because 
the mass of neutral gas in early-type gal axies varies by or- 
ders of mag nitude (e.g. Welch et al. 2010), even though their 
stellar populations are quite similar. The hybrid solution also 
removes the main weakness of the external origin hypothe- 
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sis because dust destruction is at least partially balanced by 
growth in the cold gas. As a result, the timescale that is im- 
portant to explain the demographics of dust is the cold gas 
depletion timescale, rather than the destruction timescale of 
individual grains. As the depletion timescale for the cold gas 
appears to be on the order of Gyr, largely due to star forma- 
tion, it is sufficient to maintain substantial quantites of dust in 
approximately half of all early-type galaxies. 

This hybrid solution leads to several testable predictions. 
First, the cold gas in early-type galaxies should be lower 
metallicity than the old stellar population, which may be ac- 
cessible with gas-phase abundance measurements. Second, 
the dusty early-type galaxies should have cold gas to dust ra- 
tios comparable to the accreted satellite population, unless the 
cold gas and dust are destroyed at different rates. Finally, 
there should be more evidence of mergers within the last few 
Gyr in the dusty early-type galaxies than those that appear 
dust free. 

Both external accretion and our hybrid solution also con- 
strain the lifetime of the low-luminosity AGN in the dusty 
early-type galaxies. If the only source of the cold gas and 
dust that fuels the supermassive black hole is external accre- 
tion, and these accretion events only occur every few Gyr, 
the high incidence of AGN implies that these galaxies remain 
(weakly) active for several Gyr as well, while the inactive 
galaxies remain inactive for comparable timescales. The mass 
of dust, and the two orders of magnitude more cold gas in 
the ISM, is more than sufficient to maintain accretion rates of 
< 0.01 M yr -1 for this time, as well as fuel the modest cir- 
cumnuclear star formation neccesary to produce the nuclear 
stellar disks seen in many dust- free early-type galaxies (e.g. 
iSimoes Lopes et al.ll2007f) . 

7. SUMMARY 

We have analyzed the dust in a representative sample of 38 
early-type galaxies and used these data to constrain the origin 
of their dust. This analysis is largely based on archival Spitzer 
IRAC and MIPS observations, although these galaxies also 
have HST observations that reveal the presence or absence of 
dust lanes in the central kpc, as well as high SNR, visible- 
wavelength spectroscopy that provides uniform classifications 
of nuclear activity. 

The IRAC and MIPS observations detect emission from 
dust from every galaxy that exhibits dust lanes seen in ab- 
sorption in HST observations. This dust emission is clearly 
due to both PAH particles, as detected with IRAC 8/im ob- 
servations, and cold, interstellar dust, as detected by MIPS 
observations at 70/im and 160/im. Conversely, galaxies that 
show no evidence for dust lanes in HST observations also do 
not exhibit dust emission at longer wavelengths. The absence 
of dust emission indicates that galaxies without clumpy dust 
lanes within 100s of parsecs of their centers do not have a 
substantial quantity of uniform, diffusely-distributed dust that 
would have been difficult to detect with contrast enhancement 
techniques. 

The galaxies with no evidence for interstellar dust are a 
valuable sample to estimate the stellar spectral energy distri- 
bution for old stellar populations. We use these data to esti- 
mate the ratio of the stellar emission at 8/im to that at 3.6/im 
and 4.5/im. We then use these flux ratios to scale and subtract 
the stellar contribution from the 8/im band for the galaxies 
with interstellar dust and obtain a good estimate of the emis- 
sion from PAHs in this bandpass. These flux ratios are in 
reasonable agreement with spectral synthesis models of old 



stellar populations. 

We also used the galaxies without interstellar dust to es- 
timate the ratio of the flux in the 24/im MIPS band to 
the shorter- wavelength, photospheric emission at 3.6/im and 
4.5/im. The ratio of the 24/im emission to these two shorter 
wavelength bands is fairly constant across galaxies of a wide 
range of luminosity, yet the ratio exceeds that expected from 
photospheric emission by approximately a constant factor of 
four. We conclude that the discrepancy is because emission 
from hot dust in the circumstellar envelopes of evolved stars 
makes a substantial contribution to the 24/im band. We use 
this ratio to remove the circumstellar dust contribution from 
the early-type galaxies that have dust in order to study the dust 
properties of the interstellar medium. 

A number of these galaxies have sufficient diffuse dust 
emission that we can dete rmine one or more o f the dust model 
parameters introduced by Draine &HI(5007b . These are the 
PAH fraction qpAH, the fraction of the dust heated by a strong 
radiation field 7, and the minimum starlight intensity U m i n . 
The best-fit model and the far-infrared luminosity then de- 
termine the dust mass. This analysis indicates that early-type 
galaxies are usually well fit by the same dust models that work 
well for spiral galaxies. We do find evidence that early-type 
galaxies have slightly lower values of qpAH, which may be due 
to the presence of primarily neutral PAHs and/or reflect low 
ISM metallicity. These galaxies also have lower values of 7 
and higher values of U m i n than typical spirals. 

The dust masses of these early-type galaxies are at least sev- 
eral orders of magnitude lower than the dust masses of spirals. 
The dust mass of a typical early-type galaxy detected in the 
far-infrared is 10 5-6 5 M ; however, many galaxies are unde- 
tected in the long-wavelength MIPS bands and the upper lim- 
its on the dust masses in these galaxies range up to an order of 
magnitude lower, or 10 4-5 . As our sample was designed to be 
representative of the entire early-type galaxy population, our 
results indicate that approximately 60% of typical early-type 
galaxies have > 10 5 M of interstellar dust. 

There is no correlation between the presence of dust and 
galaxy morphology. Our sample contains approximately 
equal numbers of ellipticals and lenticular galaxies and the 
same fraction of each morphological type has interstellar dust. 
We also see no strong evidence for a correlation between the 
presence of dust and environment. There is a strong correla- 
tion between the detection of dust and the presence of visible- 
wavelength emission lines, which are usually consistent with 
some form of nuclear activity. We detect dust emission from 
all of the galaxies with nuclear activity, while only approx- 
imately 25% of the galaxies classified as inactive have de- 
tectable dust emission. Based on the relatively weak nuclear 
emission in these galaxies, we conclude that the dust, and the 
larger reservoir of cold gas associated with it, is a necessary 
but not sufficient condition for accretion onto the center, su- 
permassive black hole. 

We use our dust mass measurements and the demograph- 
ics of dust in early-type galaxies to demonstrate that most of 
the dust mass can not originate from evolved stars. The main 
argument is that there are orders of magnitude variations in 
the dust mass between galaxies with the same morphological 
type and stellar mass. We also use these data to rule out an 
exclusively external origin for the dust. An external origin re- 
quires that the product of the merger-rate of gas-rich galaxies 
and the dust destruction timescale equals the fraction of dusty 
early-type galaxies. While many observations of dust lane 
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morphology and gas kinematics provide strong evidence for 
an external origin, the observed dust masses require mergers 
of gas-rich galaxies that are at least an order of magnitude too 
rare to explain the presence of dust in the majority of early- 
type galaxies. The dust destruction timescale would need to 
be nearly an order of magnitude greater than the Milky Way 
value of 4 x 10 8 years, rather than the expectation that it is 
10 4-5 years due to sputtering in hot gas, in order for purely 
external accretion to be viable. 

We propose instead that the external accretion of cold gas 
provides an environment for the continual growth of dust for 
much longer than the nominal destruction timescale of indi- 
vidual grains. If dust grains continue to grow in mass in the 
externally-accreted cold medium, the reduction in the steady- 
state dust mass should scale with the much longer depletion 
timescale of the cold gas, rather than the typical destruction 
timescale for individual grains. In fact, continual dust grain 
growth in this cold gas seems inevitable, provided the con- 
ditions are similar to those that lead to dust grain growth in 
the Milky Way. The longer gas depletion timescale of sev- 



eral gigayears is consistent with the merger rates of gas-rich 
satellites, the morphology and kinematics of the cold dust and 
gas, and the distribution of dust masses. In this scenario, the 
lifetime for the low-luminosity AGN present in most dusty, 
early-type galaxies is also comparable to this several gigayear 
timescale. 
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